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We have remeasured the absolute 1S-2S transition frequency in atomic hydrogen. A comparison with the 
result of the previous measurement performed in 1999 sets a limit of (—29 ± 57) Hz for the drift of vh with 
respect to the ground state hyperfine splitting vcs in 133 Cs. Combining this result with the recently published 
optical transition frequency in 199 Hg + against uq b and a microwave 87 Rb and 133 Cs clock comparison, we 
deduce separate limits on a/a = (—0.9 ± 2.9) x 10~ 15 yr -1 and the fractional time variation of the ratio 
of Rb and Cs nuclear magnetic moments /iRb/^Cs equal to (—0.5 ± 1.7) x 10 -15 yr" 1 . The latter provides 
information on the temporal behavior of the constant of strong interaction. 
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In the era of a rapid development of precision experimen- 
tal methods, the stability of fundamental constants becomes 
a question of basic interest. Any drift of non-gravitational 
constants is forbidden in all metric theories of gravity includ- 
ing general relativity. The basis of these theories is Einstein's 
Equivalence Principle (EEP) which states that weight is pro- 
portional to mass, and that in any local freely falling reference 
frame, the result of any non-gravitational experiment must be 
independent of time and space. This hypothesis can be proven 
only experimentally as no theory predicting the values of fun- 
damental constants exists. In contrast to metric theories, string 
theory models aiming to unify quantum mechanics and gravi- 
tation allow for, or even predict, violations of EEP. Limits on 
the variation of fundamental constants might therefore pro- 
vide important constraints on these new theoretical models. 

A recent analysis of quasar absorption spectra with red- 
shifted UV transition lines indicates a variation of the fine 
structure constant a = e 2 /Aireohc on the level of Aa/a — 
(-0.54 ± 0.12) x 10~ 5 for a redshift range (0.2 < z < 3.7) 
111. On geological timescales, a limit for the drift of a 
has been deduced from isotope abundance ratios in the nat- 
ural fission reactor of Oklo, Gabon, which operated about 
2 Gyr ago. Modeling the processes which have changed the 
isotope ratios of heavy elements gives a limit of Aa/a — 
(-0.36 ± 1.44) x 10~ 8 In these measurements, the high 
sensitivity to the time variation of a is achieved through very 
long observation times at moderate resolution for Aa. There- 
fore, they are vulnerable to systematic effects |3]. 

Laboratory experiments can reach a 10~ 15 accuracy within 
years with better controlled systematics. This type of experi- 
ment is typically based on repeated absolute frequency mea- 
surements, i.e. comparison of a transition frequency with the 
reference frequency of the ground state hyperfine transition in 



133 Cs. For an optical transition, the theoretical expression for 
the drift of its absolute frequency inevitably involves nos / I^b 
in addition to a, where //Cs is the magnetic moment of the Cs 
nucleus and /is is the Bohr magneton |4]. The microwave Rb 
and Cs clock comparison |5] even involves two nuclear mo- 
ments. The magnitude of nuclear moments results from both 
the electromagnetic and the strong interaction. 

Contributions from weak, electromagnetic, and strong in- 
teractions can be disentangled by combining several fre- 
quency measurements possessing a different sensitivity to the 
fundamental constants. In this letter, we deduce separate 
stringent limits for the drifts of the fine structure constant a, 
MCs/mb an d MRb/Vcs from combining the drifts of two op- 
tical frequencies in hydrogen and in the mercury ion with re- 
spect to the ground state hyperfine splitting in 133 Cs and the 
result of a microwave clock comparison 1 5]. Comparing mea- 
surements performed at different places and at different times 
we have to use the Lorentz and position invariance and we 
have to assume that the constants change linearly and do not 
oscillate on a year scale. With exception of this, our results 
are independent of further assumptions about a particular drift 
model and possible correlated drifts of different constants 01 . 

The experiments on the drift of the 199 Hg + 
5d 10 6s 2 S 1/2 (F = 0) -> 5d 9 6s 22 D 5/2 (F' = 2) elec- 
tric quadrupole ("clock") transition frequency j/H g were 
performed by the group of J. Bergquist at NIST between July 
2000 and December 2002 |7]. These repeated measurements 
of isng limit the fractional time variation of the ratio t'Cs/^Hg 
to (0.2±7)x lO-^yr" 1 . 

In 1999 1 8] and 2003, we have phase coherently compared 
the frequency of the (15, F = l,m F = ±1) -> (25, F' 
I . m' F =±l) two-photon transition in atomic hydrogen to the 
frequency of the ground state hyperfine splitting of 133 Cs us- 
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FIG. 1 : Simplified experimental setup for comparison of the hydro- 
gen 1S-2S transition frequency with a primary frequency standard 
(HV: high vacuum, UHV: ultra-high vacuum). A mode-locked fem- 
tosecond laser with repetition rate / rep emits the frequency comb 
which is broadened in a photonic crystal fiber, fo is the carrier- 
envelope offset frequency. 



ing a frequency comb technique |9]. The 1999 setup of the hy- 
drogen spectrometer has been described previously in iToll . so 
we present only a brief description with emphasis on the im- 
provements since 1999. A sketch of the actual setup is shown 
in FigQ 

A cw dye laser emitting near 486 nm is locked to an exter- 
nal reference cavity. The new cavity for the 2003 measure- 
ment is more stable and drifts less than 0.5 Hz/s at 486 nm. 
The linewidth of the dye laser stabilized on the new cavity has 
been characterized as 60 Hz for an averaging time of 1 s. A 
small part of the laser light is transferred to a neighboring lab 
where its absolute frequency is measured. The main part is 
frequency doubled, and the resulting radiation near 243 nm 
(corresponding to half of the 15-25 transition frequency) is 
coupled into a linear enhancement cavity inside the vacuum 
chamber of the hydrogen spectrometer. 

Hydrogen atoms from a radio-frequency (if) gas discharge 
are cooled to 5-6 K by collisions with the walls of a cop- 
per nozzle. The nozzle forms a beam of cold atomic hy- 
drogen which leaves the nozzle collinearly with the cavity 
axis and enters the interaction region between the nozzle and 
the detector shielded from stray electric fields by a Faraday 
cage. Some of the atoms are excited from the ground state 
to the metastable 25 state by Doppler-free absorption of two 
counter-propagating photons from the laser field in the en- 



hancement cavity. The 1999 measurement was performed at a 
background gas pressure of around 10~ 6 mbar in the interac- 
tion region. In the meantime, we upgraded the vacuum system 
with a differential pumping configuration. This allows us to 
vary the background gas pressure in a range of 10~ 8 -10 -7 
mbar and to reduce the background gas pressure shift and the 
corresponding uncertainty to 2 Hz. 

Due to small apertures, only atoms flying close to the cav- 
ity axis can enter the detection region where the 25 atoms 
are quenched in a small electric field and emit L Q -photons. 
The excitation light and the hydrogen beam are periodically 
blocked by two mutually phase locked choppers, and the La- 
photons are counted time-resolved only in the dark part of a 
cycle. The delay r between blocking the 243 nm radiation 
and the start of photon counting sets the upper limit for the 
atomic velocity of v < I /r, where I is the distance between 
nozzle and detector. With the help of a multi-channel scaler, 
we count all photons and sort them into 12 adjacent time bins. 
From each scan of the laser frequency over the hydrogen 15- 
25 resonance we therefore get up to 12 spectra measured with 
different delays. To correct for the second order Doppler shift, 
we use an elaborate theoretical model to fit all the delayed 
spectra of one scan simultaneously with a set of 7 fit param- 
eters ficft . The result of the fitting procedure is the 15-25 
transition frequency for a hydrogen atom at rest. 

The transition frequency depends linearly on the excitation 
light intensity due to the dynamic AC-Stark shift. We vary 
the intensity and extrapolate the transition frequency to zero 
intensity to correct for this effect ||| . 

For an absolute measurement of the 15-25 transition fre- 
quency, the dye laser frequency near 616.5 THz is phase co- 
herently compared with an atomic cesium fountain clock. We 
take advantage of the frequency comb technique recently de- 
veloped in Garching and Boulder |9] to bridge the large gap 
between the optical and radio frequency domains. A mode 
locked femtosecond (fs) laser emits a pulse train which equals 
a comb of laser modes in the frequency domain. The mode 
spacing corresponds to the repetition rate frequency / rep of 
the fs laser, and the frequency of each mode can be written 
as /„ = n/ re p + fo, where n is a large integer number and 
fo < /rep is the comb offset frequency. The repetition rate 
/ lcp can easily be detected with a photodiode. For combs 
spanning more than one octave, fo can be determined by fre- 
quency doubling an infrared mode /„ in a nonlinear crystal to 
2/„ and measuring the beat note with the appropriate mode 
fin of the blue part of the comb, yielding fo = 2f n — fa n . 

In our experiment, the frequency comb from a mode-locked 
fs Ti:Sapphire laser with a repetition rate / rep of 800 MHz is 
broadened in a photonic crystal fiber to cover more than an 
optical octave. We phase lock both degrees of freedom of the 
comb (/ rep , fo) to the radio frequency output of the Cs foun- 
tain clock to get a frequency comb with optical modes whose 
frequencies are known with the accuracy of the primary fre- 
quency standard. Any of these modes is available for optical 
measurements. 

For both the 1999 and 2003 measurements, the trans- 
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portable Cs fountain clock FOM has been installed at MPQ. 
Its stability is 1.8 x 10~ 13 r av ;!/ 2 for an averaging time r avc 
and its accuracy has been evaluated to 8 x 10~ 16 II ill at BNM- 
SYRTE. During the experiments in Garching, only a verifica- 
tion at the level of 10~ 15 has been performed. Consequently 
we attribute a conservative FOM accuracy of 2 x 10 -15 for 
these measurements. 

We have measured the 15-25-transition in atomic hydro- 
gen during 10 days in 1999 and during 12 days in 2003. For 
comparability, both data sets have been analyzed using the 
same theoretical line-shape model 1 10]. In Fig|2] the results 
of the extrapolation to zero excitation light intensity and the 
respective statistical error bars for each day are presented. The 
statistical uncertainty was significantly reduced compared to 
the 1999 measurements due to the narrower laser linewidth 
and a better signal-to-noise ratio, but the scatter of the day av- 
erages did not reduce accordingly. We ascribe this effect to a 
residual uncompensated first-order Doppler shift. Later mea- 
surements performed without the fountain clock with a delib- 
erately introduced asymmetry in the 243 nm cavity indicate an 
adjustment-dependent frequency shift (the elimination of this 
additional systematic shift should be of high concern in future 
measurements). The effect should average out after multiple 
re-adjustments of the spectrometer which have been typically 
performed twice a day. It is impossible to correct the data a 
posteriori because such details of the spectrometer adjustment 
were not recorded during the phase-coherent measurement. 
Other effects which can cause the systematic shift (intra-beam 
pressure shift, background gas pressure shift, Stark shift of the 
hyperfine levels induced by the rf gas discharge, stray electric 
fields) have been checked and can be excluded on the level of 
the observed scatter. 

The 1999 and 2003 day-dependent data were averaged 
without weighting fT3l . After accounting for the total sys- 
tematic uncertainties of 28 Hz (1999) and 23 Hz (2003) in 
the mutually equivalent evaluation processes of both measure- 
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FIG. 2: Experimental results and averages for the 1999 and 2003 
measurements of the absolute (15, F = l, rriF = ±1) — > (2S, F' = 
l,m F — ±1) transition frequency in atomic hydrogen. The back- 
ground pressure shift of +10 Hz as explained in Jg/ is added to the 
1999 value. 



ments we deduce a difference equal to (—29 ± 57) Hz over a 
time interval of 44 months. This is equivalent to a fractional 
time variation of the ratio ^cs/^h equal to (3.2±6.3) x 10~ 15 
yr- 1 . 

The frequency of any optical transition can be written as 
v = ARyF re i(a) 14 M41I . where Ry is the Rydberg energy 
and the relativistic correction F I0 \(a) takes into account rel- 
ativistic and many-body effects. F rc \ depends on the transi- 
tion in the system considered and embodies the dependence 
on a, while the parameter A is independent of fundamental 
constants. For absolute frequency measurements Ry always 
cancels out. 

Numerical calculations of the dependence of F ie \ t n g (a) for 
i/Hg on the fine structure constant a yield 11411 
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The corresponding expression for v-r gives a very weak de- 
pendence on a: 
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a— lnF ro i. H (a) ~ . 
Oa 



(2) 



The frequency of the ground state hyperfine transition in 133 Cs 
is given by 



2/Cs = A' Rya 2 — F Ie i c s (a)- 
Following fi^ll . the relativistic correction F xc \^(a) is 



a— lnF re i,cs(aO 



-0.8. 



(3) 



(4) 



Therefore, the comparison of the clock transition in Hg + 
against a primary frequency standard tests a drift of |7] 



2.0 + 0.8 + 3.2) In a , (5) 



dt z^Hg dt \ pb 
whereas the 15-25 experiment tests the drift of 



|ln^ = | fln^ + (2.0 + 0.8) ha) . (6) 
dt z/h at V Ms / 

For clarity, we set x = In a, y = In and get as 
experimental constraints 

y + 6x = (0.2 ± 7) x 10~ 15 yr" 1 (Hg+), (7a) 
y + 2.8x = (3.2 + 6.3) x 10~ 15 yr" 1 (H). (7b) 

Solving for x and y yields the separate restrictions for the 
drifts of a and pcs/l^B without any assumptions of conceiv- 
able mutual correlations (see e.g. Ja]). In this sense, this eval- 
uation is model-independent. Figure represents both equa- 
tions and their solution graphically. 

For the relative drift of the fine structure constant at the end 
of the second millennium, we deduce the limit of 

d 

x = — In a = (-0.9 ± 2.9) x 10" 15 yr" 1 . (8) 
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FIG. 3 : Possible values for the relative drift rates x and y and their 
experimental restrictions (dashed lines) as stated by eqns. (7). The 
probability density for x and y to be the true values is given by 
P(Ax, Ay) oc exp[— R(Ax, Ay)/2\ with the distances from the 
crossing point Ax and Ay and R(Ax, Ay) = (Ay+6Ax) 2 /a^ lg + 
(Ay + 2.&Ax) 2 /a^. The resulting elliptical uncertainty region, 
defined by R(Ax, Ay) = 1, gives the standard deviation when 
projected on either axis by integration over the other. The uncer- 



tainties of the measured values are an = 6.3 X 10 



yr 



and 



<j Hg = 7 X 10" 



yr 



The limit on the relative drift of /zcs/Mb is 



y 



— In — 

dt /xb 



= (0.6 ± 1.3) x 1CT 14 yr" 1 . 



(9) 



We deduce the uncertainties in expressions © and (|9j as 
projections of the ellipse (Fig. [3} on the corresponding axes. 
This is equivalent to performing Gaussian propagation of un- 
certainties when resolving Eqs. (7) for x and y. Here, the mea- 
surement results for Hg + and H are treated as uncorrelated 
even though the drift rates x and y may be correlated 0. The 
given ler-uncertainties for x and y in Eqs. (jSJl and (|9jl incor- 
porate both the statistical and systematic uncertainties of the 
hydrogen and the mercury measurements. Both limits (|SJl and 
(|9j are consistent with zero. Meanwhile E. Peik and cowork- 
ers have added a precise drift measurement on a single trapped 
Yb ion to push the overall limit even further fl3ll . 

These results allow us to deduce a restriction for the rela- 
tive drift of the nuclear magnetic moments in 87 Rb and 133 Cs. 
From 1998 to 2003, the drift of the ratio of the ground state 
hyperfine frequencies in 87 Rb and 133 Cs has been measured 
to be B 



dt u Cs 



(0.2 ±7.0) x 10~ 16 y^ 1 . 



(10) 



Substituting the corresponding dependencies F le \(a) 
these transitions yields 



for 



— In — - — 

dt v Cs dt 

Combining (|SJl, 1 11 Oi l, and 
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d_ 

dt 



Z. l n E3t = (_o.5 ± 1.7) x 10- 15 yr" 1 , 



(ID 



(12) 



where the same procedure as in Fig|5]was used with a diagram 
covering x and z = ^ ln(/XR b //xc s )- 

In conclusion, we have determined separate limits for the 
drift of a, hq s / pig and /i^b / MCs from laboratory experiments 
without assumptions of conceivable correlations among them. 
All these limits are consistent with zero. Quasar absorption 
spectra measured with the Keck/HIRES spectrograph show a 
significant deviation between the values of a today and 10 Gyr 
ago . A corresponding linear drift of a is smaller than the 
uncertainty of our result, therefore it cannot be excluded. 
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